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Abstract 



Using data collected with the "ISTRA+" spectrometer at U70 proton 
synchrotron of IHEP, we report the first measurement of the destructive 
interference in the radiative kaon decay /^~^7- We find the difference of 

the vector and axial form factors Fv-Fa = 0.126±0.027(stat) ±0.043(syst). 
The measured value is two standard deviations above the O(p^) ChPT prediction 
equal to 0.055. Inclusion of exotic tensor interaction gives Fy — Fa = 0.144 ± 
0.044(stat)±0.035(syst) and Ft = -0.0079±0.0113(stat) ±0.0073(syst), i.e. 
—0.03 < Ft < 0.01 at 90% CL, consistent both with zero and with recent 
theoretical prediction equal to \Ft\ = 0.022. 



1 Introduction 

The decay A* ^^7 proceeds via two distinct mechanisms: the internal 

Bremsstrahlung (IB) with the photon emitted by the kaon or the muon, and 
the structure-dependent (SD) decay involving emission of the photon from 
intermediate states. SD is sensitive to the electroweak structure of the kaon 
and allows for good test of theories describing hadron interactions and decays, 
like Chiral Perturbation Theory (ChPT)y4i2j. This decay is also a good place 
for searches for possible tensor interactions, see theoretical [21 IH El El E] and 
experimental jBl IHl UHl E] papers, devoted mainly to the vr — )■ ei"y decay. 
The differential probability of the decay can be written in terms of x = |^ 

and y = (where Mk is the kaon mass and and E^ are the photon 
and muon energies in the kaon rest frame). In the most general case, which 
includes hypothetical tensor term it reads: 

^ = AiBfiB{x,y) + Asd[{Fv + FAffsD+{x,y) + {Fy - Fa? fsD-{x,y)] 

-Aint[{Fv + FA)fiNT+{x,y) + {Fy - FA)fiNT-{x,y)] + ATF^fT{x,y) 
-AiBTFTfiBT{x,y) - AsdtFt{Fv - FA)fsDT{x,y). 

n , N r 1— y + r / N/ N 2xr(l — r) , , , 

fiBix.y) = ^ , r + 2 1 - X 1 - r - ^ ' ], 1 

x^(x + y — 1 — r) X + y — I — r 

fsD+ix,y) = [x + y-l -r][{x + y- 1)(1 - x) -r], (2) 
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fsD' {x,y) = [l-y + r][{l-x){l-y)+r], (3) 
fiNMx, y) = [ ^^J~J_\"_^^ m -x){l-x-y)+r], (4) 

fiNT-{x,y) = [— — r][x^ - {1 - x){l - X - y) - r], (5) 

x[x + y — 1 — r) 

fTix,y) = (x + y - 1 - r){l + r - y) (6) 

x + y-l-r rx . . 

fiBT = 1 + r ) (7) 

X X + y — I — r 

fsDT = x{l + r-y) (8) 

, ,„ V 

^SD - i- k^2q n ^2 ' 



where r = (t?^)^ with being the muon mass and 



AiB = 4r(|^)2A5D , AjNT = 4:r{^)AsD , = 4:Asd, Asdt = ^^Asd 

and AiBT = 8v^(|^)AsD. 

In these formulas Fy and are the vector and axial form factors, Ft 
is the tensor form factor. We use the prescription |3l [16] for the tensor 
interaction, a is the fine structure constant, is the charged kaon decay 
constant (155.5 ± 0.2 ± 0.8 ± 0.2 MeV [13], and Tk^^ is the width of the 
Kfj_2 decay. As in the paper[19j a minus sign precedes the interference terms, 
thus changing the sign of the form factors. Without this change in sign the 
formulas coincide with the ones given in [3l [16]. SD"^ and SD~ refer to the 
terms with different photon polarizations and do not mutually interfere. Their 
interference with IB leads to the terms labeled INT+ and INT~. The term 
SDT refers to SD~T interference, the interference with SD"*" is possible only 
in rather exotic models [Sj [6]. The interference between the tensor and 
inner bremsstrahlung amplitudes is labeled as IBT. 

The X vs y plots for different terms are illustrated in Figs.l and 2. The 
quadrangle area in these figures is restricted by the condition Xf, < x < 
Xb + 0.1 for the y interval 0.49 — 1.0 , where Xb = 1.0 — 0.5{y + \/y'^ — 4r) 
is the border of the Dalitz plot. This area, characterized by considerable 
contribution from INT~,SD^, IBT and SDT terms, and favourable background 
conditions is used in our analysis. 

Generally form factors can depend on = [p^ — p^Y = ^k{^ " 
our analysis we assume either constant form factors or the phenomenological 
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dependence for the vector and axial form factors as in [12]: Fy{q^) = Fy{0) /(I — 
gVM^) and ^^(g^) = Fa{0)/{1 - q^/Ml) with My = 0.870 GeV and Ma = 
1.270 GeV. The theoretical situation with dependence of the form factors 
is controversial. For 0(p'^) ChPT calculations [Ij there is no dependence, 
0{p^) ChPT calculations [TTf [TS] predict nearly linear dependence, similar 
to the phenomenological one. At the same time calculations [T71 [T5] in the 
framework of the light-front quark model(LFQM) give quite different dependence 
with form factors going to zero with increase. 

We assume also the g^ independence of the tensor form factor Ft- 

The absolute value of the sum of the vector and axial form factors is known 
with high precision: \Fv + Fa\ = 0.155 ± 0.008[I1], whereas the difference 
Fy — Fa is still poorly known. The latest measurent [I9] gives for Fy — Fa 
only the 90% confidence level: -0.04 < Fy - Fa < 0.24, whereas the O(p^) 
ChPT prediction is equal to 0.055 [H |2]. 

The situation with the tensor contribution is controversial. The measurements 
have been done mainly for pion decays[Hl El HDl EJ, The observation of 
the ISTRA experiment [8] was confirmed in the preliminary result by the 
PIBETA CollaborationlH Hn] having statistically significant deviation from 
the standard model without tensor term. Latest PIBETA data[ilj have however 
eliminated the evidence for the tensor term and give the limit: —0.00052 < 
Ft < 0.0004 at the 90% CL. The searches for tensor terms are done also in 
the studies of the angular distributions in the Gamow- Teller /3 decays of the 
spin polarized atoms, see for example |15| and the references therein. These 
studies are still less restrictive than 7re27 data. 

It is of interest to note that in some theoretical models |16| Ft for kaon 
decay is larger by a factor of 1/ ta.n{6c) ~ 4.4 as compared with pion decay. 
Here 6c is the Cabibbo angle. 

2 Experimental setup and event selection 

The experiment is performed at the IHEP 70 GeV proton synchrotron U70. 
The ISTRA+ spectrometer has been described in detail in papers on Kf,3 
[201 EI], Kf,3 ^ [23] and 7i-7i°7i° decays [21]. Here we recall briefly the 
characteristics relevant to our analysis. The ISTRA+ setup is located in 
the negative unseparated secondary beam line 4A of the U70. The beam 
momentum is ~ 25 GeV/c with Ap/p ~ 1.5%. The admixture of K~ in the 
beam is ~ 3%, the beam intensity is ~ 3 • 10^ per 1.9 sec U70 spill. 
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Figure 2: The same as in Fig.l for INT"*", T, IBT and SD T terms. 
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During the physics run in November-December 2001 350 milhon trigger 
events were collected with high beam intensity. This information is complemented 
by 124 M Monte Carlo (MC) events generated using GeantS [25j for the 
dominant K~ decay modes, 100 M of them are the mixture of the dominant 
decay modes with the branchings exceeding 1 % and 24 M MC events are 
the radiative K^2 decays. 

Some information on the data processing and reconstruction procedures 
is given in |20| [2T| [22| |23| [2^ . here we briefly mention the details relevant 
for present analysis. 

The muon identification (see |22| [23] ) is based on the information from 
the electromagnetic calorimeter SPi and hadron calorimeter HC. The energy 
deposition in the SPi is required to be compatible with the MIP signal in 
order to suppress charged pions and electrons. The sum of the signals in the 
HC cells associated with charged track is required to be compatible with the 
MIP signal. The muon selection is further enhanced by the requirement that 
the ratio of the HC energy in last three layers to the total HC energy 
exceeds 5 %. The used cut values are the same as in [23]. 

Events with one reconstructed charged track and one reconstructed shower 
in the calorimeter SPi are selected. 

A set of cuts is developed to suppress various backgrounds and/or to do 
data cleaning: 

0) We select events with good charged track having two reconstucted 
[x — z and y — z) projections and the number of hits in the matrix hodoscope 
MH below 3. 

1) Events with the reconstructed vertex inside the interval 400 < z < 
1600 cm are selected. 

2) The measured missing energy E^^is = Ei,eam — — is required to 
be above zero. 

3) The events with missing momentum pointing to the SPi working 
aperture are selected in order to suppress tt^tt" background ( r > 10 cm, 
here r is the distance between the impact point of the missing momentum 
and the SPi center in the SPi transverse plane). 

4) We require also the absence of the signal above the threshold in the 
calorimeter SP2 and the guard veto system GS. 

5) In order to suppress the remaining contribution at y around 0.8 we 
use the cut cos((/)(/i7)) > —0.95 where 0(/i7) is the angle between transverse 
momenta of the muon and the photon. 

We look for the signal in the distributions over the effective mass M(/i^7z/), 
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where u four-momentum is calculated using the measured missing momentum 
and assuming m.^, = 0. Effective mass spectra for the quadrangular region in 
Figs. 1 and 2 are shown in Figs. 3,4 and 5 for y-interval 0.49—1.00 with the 
step 6y = 0.03. 

The effective mass spectra have been parametrized by the sum of the 
signal and of the background. The signal form have been found from the 
signal Monte Carlo events parametrized by the sum of two Gaussians. The 
background have been found using the histogram smoothing of the MC 
background mass spectra by the HQUAD routine from the HBOOK package |26| 
This background does not ideally describe the real data, especially at low 
effective masses. This discrepance has been taken into account by addition 
of the term 

(P(4) + P(5) * M) * exp (-P(6) * M) (10) 

to the MC background. It should be noted that the contribution of this term 
in the signal regionis rather small. 

First parameter of the fit gives the number of events in the kaon peak, 
second - the position of the peak, third - normalization of the MC background. 

At small y the signal is rather small, at large y, especially in the IB 
region, it dominates over the background. The peak at the effective mass 
0.43-0.45 GeV, seen in the histogram o) is the reflection of the remaining 
Kt,2 background. 

The resulting event distribution in the interval 0.49 < y < 1.00 have been 
parametrized by function constructed using Phase Space signal MC events, 
weighted with corresponding terms (l)-(8) calculated using "true" MC x and 
y values in the corresponding Sy intervals. The results of the fit without 
tensor component are shown in Fig. 6. Here the first parameter is Fy + Fa 
fixed at the value 0.155 taken from [11], the second parameter is Fy — Fa, 
the third parameter is the tensor form factor and the fourth parameter is 
the normalization factor. In fact, the fit results are insensitive to the value 
Fy + Fa since the SD"*" and INT"*" contributions are negligible, see Fig.6. 

The fit is not perfect, but satisfactory, xVNDF= 34.11/(17 - 2) with 
around three sigma deviation from the expected and Fy — Fa = 0.126 ± 
0.027. 

Several sources of the systematic uncertainty have been studied. The 
uncertainty given by poor knowledge of the background shape was obtained 
in two ways. First, by rescaling of the errors in the effective mass distributions 
in order to have equal to one in each bin and refitting then. Second, by 
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Figure 3: Effective mass m(yU wy) spectra for the y-interval 0.49—0.67 witfi 
the step 6y = 0.03. 
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Figure 4: Effective mass m(/i wy) spectra for the y-interval 0.67—0.85 witfi 
the step 6y = 0.03. 
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Figure 5: Effective mass m(yU 2/7) spectra for tlie y-interval 0.85-1.00 witfi 
the step 6y = 0.03. 
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Figure 6: Results of fit of the event distribution with Ft = 0. 
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using polynomial parametrization instead of the form (10). First method 
leads to the deviation in the Fy — Fa equal to 0.0104, second — to the 
0.0086, the maximum of these two values was used as the contribution in the 
systematic uncertainty. The possible contribution from the z vertex position 
was obtained using increased z interval with maximum z equal to 1850 cm. 
This uncertainty is equal to 0.0329. The variation in the muon selection 
criteria gives the value equal to 0.0167. And the fit in the different y intervals 
gives the deviation equal to 0.0191. Adding the individual errors in quadrature 
we find a total systematic error of 0.0428. 

We have tried also the fit with non-zero tensor form factor. The x^/NDF= 
31.97/(17-3), Fv-Fa^ 0.144 ±0.044 ±0.035 and Ft = -0.0079 ± 0.011 ± 
0.007. 

The fit with constant vector and axial form factors gives: x^/NDF= 
33.4/(17 - 2), Fv-Fa = 0.146 ± 0.028 ± 0.046 with Ft = and xVNDF= 
31.29/(17-3), Fv-Fa = 0.136±0.060±0.045 and Ft = 0.006±0.035±0.007. 

3 Conclusions 

Our conclusions are as follows: The study of the radiative kaon decay K~2^ in 
a new region where SD~ and 1NT~ terms have maximum gives the value Fy — 
Fa = 0.126 ± 0.027(stat) ±0.043(syst). This value is 1.4 standard deviations 
above 0(p^) ChPT prediction, equal to 0.055, and probably indicates the 
need for higher order calculations or for more elaborate analysis of the 
dependence of the form factors. 

The inclusion of the tensor component gives: Fy — Fa = 0.144 ± 0.044 ± 
0.035. and Ft = -0.0079 ±0.011 ±0.007, i.e. -0.03 < Ft < 0.01 at 90% CL. 

The work is supported in part by the RFBR grant N07-02-00957(IHEP 
group). 
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